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Comment on "Viable singularity- free f(R) grav- 
ity without a cosmological constant" 

A modified f(R) gravity model has been recently pro- 
posed in [l| whose cosmological behaviour is clearly dis- 
tinguishable from ACDM. Contrary to previous opinions 
which consider that self-consistent f(R) gravity models 
distinct from ACDM are almost ruled out, the authors 
claim that the proposed model is cosmologically viable. 
Here we show that although the model satisfies some 
consistency conditions, precisely because of its departure 
from ACDM behaviour, it does not satisfy local gravity 
constraints and, in addition, the predicted matter power 
spectrum conflicts with SDSS data. 

Out of the four viability conditions imposed on f(R) 
theories 0], the proposed model satifies three of them. 
The fourth condition, namely, \fn — 1| <C 1 at recent 
epochs, is imposed by local gravity tests. Although it is 
still not clear what is the actual limit on this parame- 
ter, certain estimates @ give \Jr — 1| < 10~ 6 today, [ij. 
This condition also ensures that the cosmological evo- 
lution at late times resembles that of ACDM. However, 
in the proposed model, \fn — 1| ~ 0.2 today for a = 2 
and qo ~ —0.25. In principle, if we are only interested 
in large scales, we could ignore local gravity inconsisten- 
cies, but still the deviations from ACDM can have dras- 
tic cosmological consequences on the evolution of density 
perturbations, as discussed by several authors [1, IS0]- 

Thus, the linear evolution of matter density perturba- 
tions for sub-Hubble (krj 3> 1) modes in ACDM is given 


by the well-known expression: 
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where S = Sp/po, TL — a' /a and prime denotes deriva- 
tive with respect to conformal time r\. Notice that the 
evolution of the Fourier modes does not depend on k. 
This means that once the density contrast starts growing 
after matter-radiation equality, the mode evolution only 
changes the overall normalization of the matter power- 
spectrum P(k), but not its shape. However in f(R) the- 
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Figure 1: Linear matter power-spectra for ACDM and f(R) in 
[l[ with a = 2. Data from SDSS @. 

ories, the corresponding equation reads Q: 
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where «i = W /H 2 and n 2 = H"/H 3 . Notice the 
k 8 dependence in the last term which appears due to 
the fact that Jrr ^ 0. This means that the matter 
power-spectrum is further processed after equality and 
the transfer function is modified with respect to that of 
ACDM. This drastically changes the shape of P(k), as 
shown in Fig. 1, where normalization to WMAP3 has 
been imposed. In the figure, SDSS data from luminous 
red galaxies Q and the ACDM power spectrum from 
the linear perturbation theory with WMAP3 cosmolog- 
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ical data are also shown. Notice that ACDM gives an 
excellent fit to data with x 2 — 11 -2, whereas for the f(R) 
theory \ 2 — 178.9, i.e. 13a out. Even if we drastically 
reduced the overall normalization in a 20%, the discrep- 
ancy would still remain at the 7a level. Actually, the best 
fit would require a 32% normalization reduction and still 
would be 4.8a away (see Fig. 1). 
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